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Abstract. In the reaction of the antiproton-deuteron electromagnetic annihi- 
lation p + d — > ne + e~ the value of the invariant mass M e + e - of the e + e~ 
pair can be near or below the pp mass even for enough high momentum of inci- 
dent antiproton. This allows to access the proton electromagnetic form factors 
in the time-like region of q 2 near the pp threshold. We estimate the cross sec- 
tion da(p + d — > e + e~ n) / dM e+ e - for an antiproton beam momentum of 1.5 
GeV/c. We find that for values of M e + e - near the pp threshold this cross sec- 
tion is about 1 pb/MeV. The case of heavy nuclei 12 C, 56 Fe and 197 Au is also 
estimated. Elements of experimental feasibility are studied in the context of 
the PANDA project HE). We conclude that this process has a chance to be 
measurable at PANDA. 



1 Introduction 

The electromagnetic form factors of the proton and the neutron are basic ob- 
servables, which are the goal of extensive measurements. In the spacelike re- 
gion, /.e.for a virtual photon four-momentum squared q 2 < 0, these form factors 
give information about the spatial distribution of electric charge and magnetiza- 
tion inside the nucleon. In the timelike region (q 2 > 0) they tell us about the 
dynamics of the nucleon-antinucleon (NN) interaction. 

The theoretical models, generally based on dispersion relations G]-|5) or 
semi-phenomenological approaches [6 7 1, predict a smooth behavior of the form 
factor in the measured regions, but a peaked behavior in the timelike region be- 
low the NN threshold (0 < q 2 < Am 2 , m is the nucleon mass), due to poles in 
the amplitude (see e.g.fig.Q] taken from [ 8 1). These poles are phenomenological 
inputs, built from meson exchange, and their properties are fitted to the data in 
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the measured regions. The corresponding irregularities in form factors are re- 
lated to the transition of pp to vector mesons which can decay in e + e~ pair via 
a virtual photon. 

The mesons with a mass near the pp mass can have a quasinuclear nature, 
i.e., they can be formed by bound states and resonances in the pp system. Such 
vector mesons were predicted in the papers |9, 10 1. Note that such mesons can 
be formed not only in the pp system but in NN in general and they can have 
not only vector quantum numbers. A review on quasinuclear mesons in the NN 
system is given in ifTTl . 
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Figure 1 . Experimental data and predictions for the magnetic proton form factor in the 
domain -10 GeV 2 < q 2 < 10 GeV 2 . The figure is taken from (8). 

The under-threshold region (0 < q 2 < 4m 2 ) is called unphysical because 
it cannot be accessed experimentally by an on-shell process. Some experiments 
have been performed in the vicinity of the NN threshold, either in pp — > e + e~ 
at LEAR [12) or in the inverse channel e + e~ — > pp at Babar [13|, but they 
cannot go below this physical threshold. However, a nucleus provides nucleons 
with various momenta, in modulus and direction, and also various degrees of off- 
shellness. Therefore it offers the possibility to produce an NN electromagnetic 
annihilation with an invariant mass squared q 2 = Sp p smaller than 4m 2 . 

This possibility, which may give access to the proton form factors in the 
underthreshold region, for an off-shell nucleon, was explored in our paper ifPfl . 
The present talk is based on this paper. 

The idea to use a nucleus for that purpose was explored in the 80's using 
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deuterium 1 15 1. The reaction is then: 

pd — > e + e~n (1) 

(a crossed-channel of deuteron electrodisintegration). The aim of the present 
paper is to revive this study in view of the future antiproton facility FAIR at 
GSI. 

Other channels can give access to the off-shell nucleon form factors in the 
timelike region, including the underthreshold region; such processes have been 
studied theoretically in ref. |[T6l (jp —> pe + e~) and in refs. |fT7l[T8l (pp — > 
7r°e+e-). 

This paper is organized as follows: a theoretical study is presented in sect. [2] 
experimental aspects are discussed in sect.[3]and a conclusion is given in sect. [4] 

2 Theoretical study 

In elastic electron scattering from the nucleon e~iV — > e~N the momentum 
transfer squared q 2 = (k — k') 2 is always negative. This allows to measure the 
nucleon form factors in the space-like domain of q 2 . 

On the contrary, in the annihilation NN — > 7* — > e + e~ the mass of virtual 
photon is equal to the total cm. NN energy. Its four-momentum squared is 
always greater than 4m 2 . This allows to measure the nucleon form factors in the 
time-like domain of q 2 , above the NN threshold. In this reaction, in order to 
study the form factor behavior in a narrow domain near threshold, where non- 
trivial structures are predicted 1111 . one should have a beam of almost stopped 
antiprotons. This non-easy technical problem was solved at LEAR [12]. How- 
ever, the under-threshold domain < q 2 < 4m 2 remains kinematically unreach- 
able in this type of experiments. 
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Figure 2. The process pA — > (A — 1)7* (followed by 7* — > e + e ). 

One can penetrate in this domain of q 2 in the p annihilation on nuclei 

pA ->• (A - 1) e+e - , 
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see fig. |2 The symbol (A — 1) means not necessarily a nucleus but any system 
with the baryon number A — 1. Since extra energy of the antiproton can be ab- 
sorbed by the (A — 1) system, the e + e _ pair may be emitted with very small 
invariant mass. Therefore the two-body reaction pA — > (A — is kinemat- 
ically allowed for a very wide domain of invariant mass of the 7*, which starts 
with two times the electron mass, namely: 



One can achieve near-threshold, under-threshold and even deep-under-threshold 
values of q 2 even for fast antiprotons. This however does not mean that this 
reaction provides us direct information about the nucleon form factors. For the 
latter, we should be sure that the observed e + e~ pair (and nothing more) was 
created in the annihilation pp — > e + e~ on the proton in the nucleus, i.e., that the 
reaction mechanism is given by the diagram of fig. [3] or by a similar diagram 
where the p can rescatter before annihilation. 



Figure 3. Amplitude of the reaction pA — > (A — 1)7* in impulse approximation. 

At the same time, since the nucleons in the nucleus are off-mass-shell, the 
form factors entering the amplitude of fig. [3] are not precisely the same as found 
in the free pp annihilation. In general, the three-leg vertex F — F(p 2 ,p 2 , q 2 ) 
depends not only on the photon virtuality q 2 , but also on the nucleon ones p\ , p\. 
In the case considered, the incident antiproton is on-energy-shell: p| = m 2 , 
however the form factors depend on the proton virtuality p 2 . How the form factor 
F(pp — m 2 ,p 2 ^ m 2 ,q 2 ) vs.q 2 differs from the free one F(p? — m 2 ,p 2 = 
m 2 , q 2 ) - this depends on the dynamics determining its behavior vs.the nucleon 
leg virtuality. The nucleon form factors with off-shell nucleons were studied in 
the papers I1191I20L Generally, we can expect that the form factor dependence 
vs.p 2 is much smoother than the q 2 dependence. The p 2 dependence can be 
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determined by the nucleon self-energy corrections (i.e., by the structure of the 
nucleon), whereas the q 2 dependence in the time-like domain is governed by the 
pp interaction. The nucleon dynamics has a much larger energy scale than the 
nuclear one. The typical off-shell variation found in the papers |19 20 1 was from 
a few to 10 percent. We do not pretend to such an accuracy here. Therefore we 
neglect this effect in our calculation. We will come to this question later. In any 
case, both domains: q 2 < Am 2 , p 2 — m 2 and q 2 < 4m 2 , p 2 < m 2 are totally 
unexplored experimentally and are interesting and intriguing. 

We emphasize that though the form factor dependence on p 2 can be weak, 
the nucleon off-mass-shell effect is very important for the kinematical possibil- 
ity to reach the near- and under-threshold domain of q 2 with fast antiprotons. 
To produce the near-threshold e + e~ pairs in annihilation of a fast p on an on- 
mass-shell proton, the antiproton should meet in the nucleus a fast proton with 
parallel momentum. The probability of that, which was estimated in the Glauber 
approach, is negligibly small relative to the results presented below. However, if 
the effective mass p 2 of the virtual proton is smaller than m 2 (that is just the case 
in a nucleus), then the near- and under-threshold e + e~ pairs can be produced in 
collisions with not so fast intra-nucleus nucleons. This effect considerably in- 
creases the cross section. To have an idea of the order of magnitude which one 
can expect for this cross section, we will calculate it in the impulse approxi- 
mation. Numerical applications will be done for the lowest antiproton beam 
momentum foreseen in future projects. Namely, at the High Energy Storage 
Ring at FAIR-GSI this value is 1.5 GeV/c. 

2.1 Cross section calculation 

At first, we consider the case of the deuteron target. If we know the amplitude 
of the reaction pd — > e + e~n: Mp d _> e + e - n (to be calculated below), then the 
corresponding cross section is given by: 



dcTpd-te+e-n ~ " |j \Mpd-^e+e-n\ 2 (2) 

(4) d 3 p e + d 3 p e - d 3 p n 

XS { >(pp+Pd-Pe+ ~Pe- -Pn) , , 3o , x 3o , Q ^ 3o 

[2n)- i 2e e + (2ir) i 2e e - (27r) d 2e„ 

where / results from the flux factors. Here and below we imply the sum over the 
final spin projections and average over the initial ones. 

Our estimations are carried out in the impulse approximation, when the 
mechanism is given by the diagram of fig. [3] Then the total amplitude squared 
\Mp d ^ e + e - n \ 2 is proportional to the annihilation amplitude squared \Mp p ^ e + e - | 



5 



V.A. Karmanov, H. Fonvieille 



and to the square of the deuteron wave function \ip\ 2 : 

\M pd ^ e+e - n \ 2 = 4m \M pp ^ e+e - | 2 |^| 2 , (3) 

and | ip\ 2 is normalized to 1. 

We are interested in the distribution in the invariant mass A4 of the final 
e + e~ system. To find it, for fixed value of A4, we can integrate, in some limits, 
over the angles of the recoil neutron (determining the neutron recoil momentum) 
and over the angles of the emitted e + e~ in their center of mass. This can be 
done using standard phase volume techniques. The derivation is presented in 
detail in 03). The final result reads: 

rf ^ +C ~" = Ofcp-e+e- (M) V(M), (4) 

where rj(Ai) is the distribution (given by eq. (O below) of the e + e~ invariant 
mass A4 and <J pp ^ e + e - is the cross section of the pp —> e + e~ annihilation 
at the total energy A4 . 

The calculation of cr pp ^ e + e - is standard. To estimate the nuclear effect, we 
omit the nucleon electromagnetic form factors. Then the pp —> e + e~ cross 
section obtains the form: 

_ 2a 2 7r(2m 2 + M 2 ) 
^ e+e ~ " 3M 2 mp pMb ■ (5) 

Though Cpp^e+e- depends on VW, the main (nuclear) effect is determined 
by the factor rj(M) lfl4l: 

r)(M) = /0 7 ? " / \m?dz. (6) 



(27r) 2 v /s~^ 



i 



Here z — cos 8, where 6 is the angle, in the cm. frame of the reaction, between 
the initial deuteron momentum p* d and the final neutron momentum p^. The 
argument of the wave function k depends on z. This explicit dependence is 
given in fl4"l . 

Since, as mentioned, the wave function squared |7/>(fc)| 2 is normalized to 1, 
the distribution r](M.) is also automatically normalized to 1: 

r](M)dM = 1. (7) 

The total cross section is obtained by integrating (O in the finite limits 

Mmin < M < M max , where M min = 2m e , M max = ifsp~d - m. Ne- 
glecting the electron mass, we can put M.min = 0. For p p j a b = 1.5 GeV/c the 
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value M- max is high enough and provides the normalization condition (|7]i with 
very high accuracy. 

To emphasize more distinctly the effect of the nuclear target, we can rep- 
resent the cross section © of the annihilation pp — > e + e~ on a free proton 
similarly to eq. (|4]l: 



where ap p ^ e + e - is defined in ©, s pp — (p p + pp) 2 . The fact that in the anni- 
hilation on a free proton the mass of the final e + e~ pair is fixed is reflected in 
dHJ in the presence of the delta-function. Comparing this formula with we 
see that the effect of the nuclear target results in a dilation of the infinitely sharp 
distribution S(Ai — yjs pp ) in a distribution of finite width rj(Ai). The dilation 
of a distribution does not change its normalization: i](Ai) remains normalized 
to 1. 

2.2 Analysis and numerical calculations 

At first glance, the small near-threshold cm. energy Ai « 2m in the collision 
of a fast p (pp = 1500 MeV/c) is achieved, when the antiproton meets in the 
deuteron a fast proton having the same momentum as the p, in modulus and 
direction. The protons with such a high momentum are very seldom in deuteron. 
For this mechanism, the cross section would be very small. However, the near- 
threshold value of Ai is obtained in other kinematics. As we mentioned, the 
proton momenta k in the deuteron wave function ijj(k) in eq. © start with 
k « k m i„ ~ 360 MeV/c only (that corresponds to z « 1). The main reason 
which allows to obtain in this collision the value Ai ~ 2m is the off-shellness 
of the proton: m* < 0.85m instead of rn* = m. This 15% decrease relative to 
the free proton mass is enough to obtain the invariant pp mass Ai ~ 2m, when 
one has the two parallel momenta: 1500 MeV/c for p and 360 MeV/c for p. 

The cross section d(Tp d ^ e + e - n /dAi, eq. ©, has been calculated for an 
antiproton of momentum p p — 1500 MeV/c on a deuteron nucleus at rest, with 
the deuteron wave function ||2~T1 . incorporating two components corresponding 
to S- and D-waves. The result is shown in fig. [4] The maximum of the cross 
section is at Ai = 2257 MeV, that corresponds to the p interacting with a proton 
at rest (and on-shell). The cross section integrated over Ai is equal to 43 nb. We 
remind that these calculations do not take into account the proton form factor. 
Its influence will be estimated below. The numerical integral over Ai of the 
function r)(Ai), eq. ©, is ps 1, in accordance with the normalization condition 




da. 



(8) 



0. 
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~H + — — -t- — 

Figure 4. The cross section — P Xvi " — ™ °f tne reac ti° n pd — > e e n VJ..M, in the 
interval: 1750 MeV < M < 2350 MeV, calculated for a pointlike proton. The arrow 
indicates the pp threshold. 



The pp threshold value M = 1880 MeV is on the tail of the distribution, 
far from the maximum. Relative to the maximum, the cross section at threshold 
decreases approximately by a factor 600. The numerical value at the threshold 

is: 

da 
dM 

In fig.|5]this cross section is shown in the near-threshold interval 1830 MeV < 
M < 1930 MeV. The integral over M in a bin of width 100 MeV centered on 
the threshold is: 

,.1930 MeV j 

aa pd ^ e+e - n dM ^ 10Q pb 



pb 

M=2 m = 1 MeV' 



' 1830 MeV dM 

These estimations take into account the suppression resulting from the mo- 
mentum distribution in deuteron. However, they do not incorporate the form 
factors of the nucleon. To incorporate them in a simplified way, one can con- 
sider an effective form factor \F\ which depends on M, and include it in the 
integral: 



Vpd^e+ e -n = J <Jp P ^e-(M) T,{M) \F{M)\'dM (10) 

where <Jp p ^, e + e - (M) is the cross section for pointlike nucleons given in eq.©. 
To have an estimate of this integral, we have taken the effective proton form 
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Figure 5. The same as in fig. [4] but near the pp threshold, in the interval 1830 MeV < 
M < 1930 MeV. 



factor measured in ref. iflD . By doing this, we neglect all off-shell effects. We 
interpolate linearly between the measured values, and we limit the inte- 

gral to the region M. > 2m. In this way we obtain o-pd^e+e-n — 1 n t>, which is 
comparable to the total cross section <7p P ^ e + e - on a free proton at A4 = 2257 
MeV. We point out that at threshold, our differential cross section dapd^ne+e- 
of 1 pb/MeV (eq.©) is not suppressed by any factor, since there the form factor 
|F| seems to be close to 1 experimentally 122II23II . Below this threshold one may 
expect a form factor effect larger than one. 

2.3 Annihilation on heavier nuclei 

For A > 2, we should take into account the possibility of excitation and breakup 
of the final nucleus A — 1 in the process pA — > (A — 1)7*. The result con- 
tains the sum over the final energies of the residual nucleus and the integral over 
a continuous spectrum. That is, the function |-0(A:)| 2 in eq. (|6]l is replaced by 
the integral /J' max S(E, k)dE, where S(E, k) is the nucleus spectral function 
giving the probability to find in the final state the nucleon with the relative mo- 
mentum k and the residual nucleus with energy E. For high incident energy we 
can replace the upper limit by infinity. Then we obtain: 

/ S(E,k)dE = n{k), 

JE min 
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where n(k) is the momentum distribution in the nucleus. 

To estimate the cross section on heavy nuclei, we will still use eqs. ©, (O 
but with the two following changes, (z) We replace the deuteron momentum 
distribution by the nuclear one. (ii) We multiply (O by the number of protons Z. 

The numerical calculations were carried out for the 12 C, 56 Fe and 197 Au nu- 
clei with the nuclear momentum distributions found in the papers by A.N. Anto- 
nov et al: |24][25]. Near threshold, z'.e.at M = 1880 MeV, for all three nuclei 
we obtain very close results given by (compare with eq. (O for deuteron): 

da pA ^ +e - x w65Z _gh_ (n) 
dM MeV 

Multiplying by the charge Z (Z( 12 C) = 6, Z( 56 Fe) = 26, Z{ 197 Au) = 79) and 
integrating (fTTT ) over a 1 MeV interval near Ai = 1880 MeV, we get: 

o-( 12 C) = 39 pb, cr( 56 Fe) = 0.17 nb, cr( 197 Au) = 0.5 nb. 

These results were obtained without taking into account the absorption of p 
in nucleus before electromagnetic annihilation. This absorption was estimated 
in Glauber approach. It reduces these cross sections by only a factor 2. 

2.4 Beyond the impulse approximation 

There exist other possible mechanisms for the process pd —> e + e~n. One of 
them is the initial state interaction, which includes rescattering (not only elastic) 
of the initial p in the target nucleus. In the rescattering, the incident p looses 
energy and therefore the proton momentum needed to form the invariant mass 
A4 R3 2m becomes smaller. The probability to find such a proton in deuteron is 
higher. Therefore initial state interaction increases the cross section. 

Other processes are discussed in [14 |. We emphasize that in any case, what- 
ever the intermediate steps are in process ([T), the e + e~ pair of the final state must 
come necessarily from the baryon-antibaryon electromagnetic annihilation, pp 
or fin, because there is only one neutron left at the end. It cannot come from 
another process, even if there are complicated intermediate steps, like rescatter- 
ing, etc. Therefore this e + e~ pair is a direct and very little distorted probe of the 
baryon-antibaryon electromagnetic annihilation vertex. 

3 Experimental aspects 

Experimental aspects have been investigated in detail in [14| in the case of a 
deuteron target, i.e. for the three-body process p d 
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the PANDA project at FAIR-GSI: an antiproton beam momentum of 1 .5 Ge V/c 
and the detection of the lepton pair. The count rate in the near-threshold region 
of AA is small but not negligible. The main difficulty is to identify the reaction 
among the hadronic background which is about six orders of magnitude higher. 
First elements of strategy were presented for this background rejection, based 
on particle identification, detector hermeticity, and missing mass resolution. 

One should also also note that the luminosity in (pA) decreases with the 
atomic charge Z of the target nucleus 0], in a way that roughly compensates the 
increase of cross section with Z reported in sect. 12.31 

Although the subject would require a much more detailed study, it was con- 
cluded that this process has a chance to be measurable in PANDA, given the 
very good design performances of the detector. 

The antiproton momentum 1 .5 Ge V/c just corresponds to the threshold value 
of creation of the AA pair on a free proton. Therefore the virtual creation of the 
AA pair in reaction ((TJ is not suppressed by the nucleon momentum distribution 
in deuteron and contributes just in the domain of the peak of fig. [4] that allows 
one to study the AA threshold region with good statistics. In the AA system, 
the quasi-nuclear states were predicted in 11261 and, similarly to the NN quasi- 
nuclear states, they should manifest themselves as irregularities in the cross sec- 
tion. The contribution of the channel pp — > AA — > e + e~ in the total cross 
section pp —> e + e~ (the latter equals 1 nb, see sect. 12.21 above) was estimated 
in 1 27 1 as 0. 1 nb, i.e. 10% of the total cross section. Therefore we expect that the 
structures caused by the channel pp — > AA — > e + e~ can be observed in process 
(1) in the region of mass AA near the A A threshold. 

4 Conclusion 

We have studied the reaction pA —> (A — 1)7* (followed by 7* — > e + e~). This 
process gives access to the pp annihilation pp — > 7* at invariant masses y/~Sp p 
which are below the physical threshold of 2m, due to the proton off-shellness in 
the nucleus. In this way a possibility exists to access the proton timelike form 
factors in the near-threshold and the totally unexplored under-threshold region, 
where NN bound states are predicted. 

The differential cross section da/dAA has been calculated as a function of 
the dilepton invariant mass AA, for an incident antiproton of 1.5 GeV/c momen- 
tum on a deuteron target (and heavier nuclei). We find that the pp threshold 
(AA = 2m) is reached for a minimal proton momentum fc m j n =360 Me V/c in the 
nucleus, and at this point the cross section is about 1 pb/MeV. 

The estimations carried out above show that this process has a chance to 
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be measurable in PANDA, provided the very good design performances of the 
detector. 
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